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1. Introduction 

The Dyson-Schwinger equations of Coulomb gauge quantum chromodynamics (QCD) rep- 
resent one of the many techniques being currently explored in the hope of one day being able to 
describe confinement and the hadron spectrum from first principles. As with most difficult prob- 
lems, it is useful to be able to compare and contrast different approaches to gain further insight. 
One aim of this talk is to compare a leading order truncation of the Dyson-Schwinger equations 
||l|] to the gap equations for the static gluon and quark propagators obtained within a quasi-particle 
approximation to the canonical Hamiltonian approach [Q, ^, ^]. 

The talk starts with a brief review of Coulomb gauge within the first order formalism, including 
a discussion of the charge constraint that emerges from the incompleteness of the gauge fixing. To 
avoid problems stemming from the nonlocality of this formalism, an Ansatz is introduced such that 
the Dyson-Schwinger equations can be derived. The reduction of the truncated Dyson-Schwinger 
equations to the gap equations for the static propagators and the link to heavy quarks will be shown. 
How the heavy quark limit provides an intuitive explanation for the charge constraint and infrared 
divergences as being unobservable constant shifts in the potential will be discussed. 



2. Coulomb gauge in the first order formalism 



Let us begin by considering the functional integral associated with QCD (in Minkowski space): 



m 



^^qp, + ^E^-E^-^B"rBl;\ (2.1) 



where generically represents the functional integration measure over all fields present. The 
(conjugate) quark field is (q) q^^ where the fundamental color, spin and flavor indices are denoted 
collectively with the index j8 and the position argument with subscript x. The Dirac y-matrices obey 
the usual Clifford algebra {y^ = 2g^^ with metric g^^ = diag(l, — 1) (we explicitly extract all 
the minus signs associated with the metric such that all components of a spatial vector x are written 
with subscripts, i.e., xi). The temporal and spatial components of the covariant derivative in the 
fundamental color representation are given by 



'Ox 



-iga^r, D, = V, + igA",T'' 



(2.2) 



where a" (= A^") and A" are the temporal and spatial components of the gluon field, respectively, 
and where the superscript a denotes the color index in the adjoint representation. The generators 
obey [7*^,7^] = if"^'^T'^, where the f'^^'^ are the structure constants and we use the normalization 
jj-j jfo-j _ §ah xhe chromoelectric and chromomagnetic fields are written in terms of the gluon 
field as 



BqA", - Dfa^, = v., X A? - ^gf'A'^ x A 



1 



■ahc ~Ah 



with the spatial component of the covariant derivative in the adjoint representation given by 



D 



ah 



5"''V,--gf"'A 



■acb 'Xc 

X' 



(2.3) 



(2.4) 



The QCD action is invariant under gauge transforms A — )• A^ = UAU^ — i/g{dU)U\ q^ q^ = 
Uq, where Ux = exp{— iO^T"} is a spacetime element of the SU (Nc) group parametrized by 6^. 
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Because of this invariance, the functional integral contains a divergence due to the integration over 
the gauge group. When calculating Green's functions such as propagators, it is thus necessary to 
fix the gauge and our choice is Coulomb gauge: V • A = 0. The Faddeev-Popov (FP) technique to 
fix the gauge involves inserting the identity 



Det 



ah I 



5F" 









(2.5) 



F=0 



into the functional integral. However, in Coulomb gauge where f = V • A and the FP kernel reads 
M{x,y) ~ — •Dx5{x — y), there is an obvious problem when the gauge transform parameter 6" 
is spatially independent: 

-V,-Dfd\xQ)=Q (2.6) 

(there are no temporal derivatives) such that the FP determinant automatically vanishes. Coulomb 
gauge is incomplete in this respect. The resolution of the temporal zero modes of the FP operator 
leads to a constraint on the total color charge of the system in the first order formalism [^. The 
identity, Eq. (2^), is modified to 

1 = j ^65 [l 



Det 



M''''{x,y) 



(2.7) 



where &6 and Det explicitly exclude the temporal zero modes, d{xo). The Coulomb gauge fixed 
functional integral is thus 



V- A Det 



-VD 



(2.8) 



The conversion to the first order formalism goes as follows [^, ^ |5|]. An auxiliary vector field (tt) 
is introduced via 



1 



exp <j I j dx^E"-E° 



1 



Sjtz exp J dx 

and split up into components {(j) is the longitudinal part of n) with 

const = J &(j)9rexpl-i J dxz^ (v^n^ + Vl(^^ 



(2.9) 



(2.10) 



Changing variables tt — tt — V0 and integrating out the Lagrange multiplier, the functional integral 
now has the form 



^<I>5(V-A)5(V-7r)Det 



-VD 



1.9" 



(2.11) 



where the action, .5/" , is at most linear in the temporal gauge field, a: the corresponding term is 

= I dxat {y.-Df^^+gf'^A^-7i',+gq„Afnam-^) ■ (2-12) 

(Incidentally, the above form of the functional integral is the starting point for studying perturbation 
theory in the first order formalism ^, ^.) Importantly, the a-field can be integrated out to give 



Z= I &<^5(v -Ajsiv- TiJ Det -V D 5(v-D<l)+p^e 



j'y") 



(2.13) 
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where the color charge, p, includes both gluonic and quark contributions: 



(2.14) 



The field can be integrated out by using the eigenfunctions of the Faddeev-Popov operator as a 
complete orthonormal basis for an expansion, the crucial point being that one must remember the 
temporal zero modes Including the -dependent part of the action, the explicit expression is 



&<^8 ( V-D(/) + p ) exp <; - / dx^y^cj): 



dxp I Det 



-VD 



exp 



where 



-Vv-d; 



cb 



(2.15) 



(2.16) 



Notice the appearance of the inverse (modified) FP determinant, that will cancel against the original 
in the functional integral. The 5-functional constraint that emerges constrains the total color charge, 
the spatial integral arising from the projection onto the temporal zero mode. In order to study its 
effect, we rewrite this 5-functional constraint in Gaussian form 



dxp 



~ lim ^('^)exp<' -- / dxdyp"{x)'^5"''5{xQ-yo)p\y) 



(2.17) 



where is a constant, .yV(^) is a normalization factor to be included implicitly in the functional 
integral measure, and the limit — )• oo will be taken only at the end of any calculation. With this, 
our functional integral now reads 



V-A 5 V-TT e 



(2.18) 



with the action 



,5^ = dx< q^^ i^/^dox + iy-Dx — m 



IPx 



- / dxdyp-J"\xj)p'y 



and where F is the Coulomb kernel, but shifted by a spatial constant proportional to 



F"\x,y) = Ff5{x-y)+^8"'8{xQ -yo). 



ab ; 



(2.19) 



(2.20) 



There exists a useful connection between the Coulomb kernel and the temporal gluon propagator 



[|10p. Redoing the analysis for the functional integral in the presence of a source (p) for the temporal 
gluon field, the temporal gluon propagator is defined as 



^'^'^^'^'^^ Z[p] 5ip^5ip,^ 



(2.21) 



p=0 



The presence of the source p only alters the above action, Eq. ( 2. 19| ), by replacing p with p 



p + p. As noted [10], the temporal gluon propagator has a purely instantaneous part given by the 
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expectation value of the Coulomb kernel since it involves only spatial derivatives. In our case, 
where the kernel is shifted by a constant, we see that 

Kai^^y) ^<iFf5{x-y)+i'^6"''> d{xo -jo) +non-inst. (2.22) 

To recap, by writing the Coulomb gauge functional integral in the first order formalism, the FP 
determinant cancels after integrating out the temporal and longitudinal fields and Coulomb gauge is 
thus ghost-free [^, ^. What remains of the gluon field are the two transverse vector components A 
and n (which would give rise to the two polarization states of photons in quantum electrodynamics). 
Treating the temporal zero modes of the FP operator explicitly, it is further seen that the total color 
charge must be conserved and vanishing [^]. This is nothing more than the application of Gauss' 
law. Writing the total charge constraint in Gaussian form, the Coulomb kernel is shifted by a spatial 
constant - eventually however, we must take the limit where this constant diverges. We shall see 
though that this is not a problem in the end. 



3. Truncated Dyson- Sch winger equations 

Having written down our functional integral in the first order formalism, we would like to use 
it. Unfortunately, the Coulomb kernel term (F) is nonlocal because of the presence of the inverse 
FP operator. In order to derive Dyson-Schwinger equations, we therefore make a truncation Ansatz 
whereby we replace the Coulomb kernel with its expectation value ^ : 

F"\x,y) ^ [F{x-y)+'^]d'"d{xo-yo) (3.1) 

where F is now some purely spatial, scalar function which will serve as nonperturbative input into 
the system. Note that this Ansatz still includes the tree-level term, such that one-loop perturbative 
results could still be obtained at this stage. The action is now local, and given the form of the color 



charge p, Eq. (2.14), the Coulomb interaction term pFp now involves a set of effective four-point 
vertices (see below for their exphcit form). In effect, by converting to the first order formalism, 
we replace the dynamics of the nonperturbative towers of Dyson-Schwinger equations [Q] and 
Slavnov-Taylor identities | [IT| ] involving the temporal (a), longitudinal (0) and ghost degrees of 
freedom with our leading order Ansatz for F. 

Since we have only modified the Coulomb interaction part of the action, many of the propa- 
gator and vertex Green's functions in the present formalism can be read off from previous studies 
[^, ^]. In particular, the propagators (W in our notation) in momentum space are given by: 



A,{k) ' 

-Y-kAjk)+Bjk) + ')^knY-kAjk) 



W^,ap{k) = -j^[fkoA,{k)-Y-kA,s{k)+B^{k) + fkofkA4k) 



Ag{k) = k'onA>^)-'i^AA{k)r„Ali) + iO+, 

Af{k) = klA^{k)-k^Al{k)-Bi{k) + klk^Al{k) + iO+, (3.2) 
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where the various dressing functions arise from the decompositions of the proper two-point func- 
tions (F): 



nnij 



(k) 



Anij 

T^ah 
A AAij 

.(_0) 
qqaji 



{k) 
(k) 



{k) + iij{k)r„„{k) 

5ijrAn{k) + hj{k)TAn{k) 

ij{k)YAA{k)+lijCk)TAA{k) 

fkoAt{k) - fkA,,ik) - B,nik) + fkofkAdik) 



(3.3) 



In the above, /,y and tij are the usual longitudinal and transverse spatial projectors, respectively. The 
components of the gluon propagator are spatially transverse because we are in Coulomb gauge. The 
dressing functions are all scalar functions of kg and Ir separately, due to the noncovariance. At tree- 



level r„ 



- An 



AA 



■A, 



1 , B„, = m and all others vanish. Notice the matrix inversion 



structure of the components of the gluonic and quark propagators, with the denominator factors 
and Af - these will turn out to be important. The tree-level quark-gluon (F^^a). three- (Taaa) and 
four-gluon (Faaaa) vertices are also unaltered from [^, ^], although their explicit form will not be 
needed here. With our Ansatz to replace the Coulomb kernel with its expectation value, the new 
tree-level vertices explicitly read (all momenta incoming) 



r'ff''5n5jkF'Hki +k,) + f''ff"'djidikF'^f{k, +ks) 



ri^L Q s:{k\,k2,kT,,kA] 



ig^r''[fTf]^pd,-jF'^f{k,+k2), 



r/3 



(3.4) 



With a little practice, the Dyson-Schwinger equations are not difficult to derive (although keep- 
ing track of the signs when quarks are present is somewhat tedious). Generically, their structure 
arises from the Legendre transform and repeated functional differentiation of the generating func- 
tional, giving the characteristic sequence of loop integrals. Such a derivation in Coulomb gauge is 
given in Refs. ||7|, |l^, ^. Omitting the two-loop contributions, the Dyson-Schwinger equations for 
the proper two-point functions, in the system considered here [p, are presented diagrammatically 
in Figs. [T] and |2[ Because of the existence of the mixed gluon propagator Wak, certain loops 
involve a sum over the two gluonic field types A, n which is denoted by B and C in the diagrams. 
In addition to the truncation to omit two-loop contributions, we further restrict to considering only 
those terms arising from the tree-level four-point vertices involving the Coulomb kernel F, i.e., 
we throw away the FqqA, ^aaa and F^^^ tree-level vertices. The remaining loops of the Dyson- 
Schwinger equations are thus tadpole contributions involving the propagators and our input Ansatz 
for F, forming a closed set of equations. The input we have in mind is motivated by the connection 
to the instantaneous part of the temporal gluon propagator. In momentum space and omitting the 
perturbative contributions, we will assume the strongly infrared enhanced form: 

d"''F{P) + 5"V(27r)^5(^), = SnOc/k^ 



g^CpF"^ 



(k) 



(3.5) 



where Cp = {N^ — \ )/2Nc and Oc is the Coulomb string tension |[T3|]. Note that g^F is a renormal- 
ization group invariant quantity in Coulomb gauge |I0|], so is ideal for use as input. We shall 
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Figure 1: Dyson-Schwinger equations for Tim, F^j^ and T^a, omitting two-loop terms. Wavy lines denote 
proper functions, the large filled blob indicating the dressed function. Springs denote gluonic propagators, 
lines denote the quark propagator and all internal propagators are dressed. Small blobs indicate tree-level 
vertices and large circles denote dressed vertices. The gluonic field types B and C denote the sum over A and 
n contributions arising due to the presence of mixed gluon propagators. See text for details. 



show that this truncation results in gap equations for the static gluon and quark propagator dressing 
functions that can be compared to those derived in the canonical Hamiltonian approach, Refs. [Q] 
and iQ], respectively. 



4. Leading order static gluon equation 

Let us consider the truncated Dyson-Schwinger equation for the mixed gluonic proper two- 
point function, Yjia (middle line of Fig. [T]). Recognizing that the color structure of the quark tadpole 
loop vanishes, the equation can be written (dco = d"^ (O / {2nY) 

8ijTAn{k) + lij{k)TAn{k) = 8ij-ig^N, I ^^^^^^^^Uj{cb)F{k-(o) (4.1) 

J l<.Ql\g\(X)) 

where we have expanded the two-point functions using Eqs. (|3.2y3.3D and the four-point function. 



Eq. (3.4), subsequently resolving the color structure. Since F is energy independent (coming from 
the instantaneous Coulomb kernel) and all dressing functions are even functions of energy, the 
energy integral of the above is overall odd and vanishes. This immediately gives the results 

TAn = l, TAn = 0, ^g{k) = kl -k^T AA{k)Y ^^{k) + lO + . (4.2) 
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Figure 2: Dyson-Schwinger equation for the quark two-point function, omitting two-loop terms. On the left- 
hand side, the filled blob indicates the dressed (inverse) propagator, otherwise notation is as in the previous 
figure. See text for details. 

Turning to the truncated Dyson-Schwinger equations for Tj^i and Taa (first and last lines of Fig. |I|, 
respectively), after sorting out the decompositions and color factors as above, we have (the longi- 
tudinal dressing functions Yjui and Yaa play no role here) 



TAA{k) 



-F{k-(o)tji{k)tij{a)), 



d (0 a)^rAA{(o] 



■F{k-(o)tji{k)tij{(d). 



(4.3) 



The energy integrals do not involve ko because of the energy independence of F and we thus see 
that r^j;!: and Taa are energy independent. The energy dependence of the denominator factor 
is now reduced such that we can now define the static (i.e., energy integrated or equaltime) gluon 
propagator in terms of a single dressing function, G [|I|] : 

''^''^Wfl:Ak) = 5"'tn(k)-L 

2\k\ 



AAij 



= / ^W^,ik) = d%{k)^G0?^'l' 



(4.4) 



A similar expression exists for the static Ti-propagator, Wji^ ■ Further inserting the Ansatz form for 



F from Eq. ( [3.5| ) and writing the spatial integrals in terms of G, the equations become 



TAAik" 



"^^F \/ k'^ ■' y 

N, 



4Cf 
N, 



2Ce- . /z,2 4Ce- J 



deb ap- 

~7^p 



G{a)^)-'/^F{k - a>)tji(k)tij{a)) (4.5) 



where deb = deb/{27l)^. The proper dressing functions Fjck and Taa have contributions linear in 
(the constant that arises from the charge conservation) and also from the potentially infrared 
divergent spatial integrals over F ~ — ebY, if we use Eq. (3^) as input. However, further 
utilizing the definition of G, Eq. (^!^, we find that we can combine the above coupled equations 
into a single equation for the static gluon propagator dressing function: 



G{k^ 



1 + 



1 

4Cf 



deb 



eo^ 



F{k-(b)tji{k)tij{(b) 



G{cb') 



2^1/2, 



P G(«2)l/2 



(4.6) 
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This is the gluon gap equation and is identical to that originally derived from the canonical approach 
The troublesome terms proportional to 'io drop out and the infrared divergence of the spatial 
integrals is canceled (this is explicitly verified in [p). It thus appears that under this (leading order) 
truncation, the static gluon propagator contains the physical dynamics of the system whereas the 
full propagator (in particular, its pole position) is unphysical. We shall discuss this at the end of 
the next section. It is known that for an interaction of the type given by Eq. (3^), the solution 
to Eq. (4^6) is of the massive type ^, in contradiction to the expected Gribov type solution 
[14]. However, from the canonical approach, it is known that the gap equation receives significant 
infrared contributions from the ghost loop ('curvature') [||] which is missing from the leading order 
truncation presented here. 



5. Leading order quarks and the heavy Umit 



The analysis for the quark Dyson-Schwinger equation is very similar to that previously de- 
scribed for the gluon. This similarity arises because the color charge p, Eq. ( |2.14 ), treats the 
gluonic and quark contributions on an equal footing. Truncating the equation (Fig. inserting 
the appropriate factors, Eqs. (3.2,3.3|3|4), resolving the color factors and projecting out the Dirac 
components, we obtain four coupled equations (one for each of the dressing functions). Two are 
trivial because they involve odd energy integrals: 



A,{k) = l-ig^Cp 



d (o (OqA,{co)F {k - (o) 
koAf{(o) 

ii(0(Ook-(bAd{co)F{k-(o) 



kok^Af{(o) 



such that 



At = l, Ad = 0, Af{k) = kl-k^A^{k)-Bl{k) + iO^ 
The other two equations are 



(5.1) 



(5.2) 



A,{k) = l + ig^Cp 



dcok-(bAs{(o)F{k- co) 

pAfico) 
dcoB,„{co)F{k- (o) 



(5.3) 



where the energy independence of F again means that A^ and B„, are purely spatial. As for the 
gluon, we can now write the static quark propagator in terms of a single dressing function, M: 



Y-k-M{P) 



aj3 



2^P+M{Py 



M{P) = 

A,.{P) 



(5.4) 
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Inserting the form of F given by Eq. ( |33| ), these equations become 



1 + 1 



2 



^ 1 /■ da)k-(oF{k- (b) 



1 "^MfP 



_^ 1 /■ d(bM{(o^)F{k-(b) 



(5.5) 



Like for the gluon, the dressing functions for the quark propagator are also linear in and in- 
volve potentially infrared divergent spatial integrals. However, the coupled equations can also be 
combined into a single gap equation for M: 



Mikn =m + - ^ 

2J y/W+MW? 



(5.6) 



where the '^-dependence and the potential infrared divergence cancels (also verified explicitly in 
Ref. [|l]]). This equation is very well-known as the Adler-Davis truncation [Q] and was originally 



derived using the canonical approach. With the interaction, Eq. (3.5), the solution does exhibit dy 



namical chiral symmetry breaking [|l], Q], but to a quantitatively too small degree for phenomenol- 
ogy - this leading order truncation requires further contributions [p^. 

It is possible to make the connection between the truncated quark propagator and the known 
Coulomb gauge heavy quark limit [ [I^ ] (see also Carina Popovici's contribution to these proceed- 



ings). This is done via a spin-decomposition of the full quark propagator []17[]. Introducing the spin 
projectors 

P± = {l±f)/2, P++P- = \ (5.7) 



the full quark propagator can be written as 



W^^a^ik) = [{P++P^)W^,{k){P+ +P-)]^p , 



(5.8) 



We now consider the heavy quark limit in the Coulomb gauge rest frame: \k\/m — 0. Using 
Eq. ( |5.6[ ), we can make an estimate for the static dressing function M (this is confirmed numerically 
in [|l|]). The function F{k — G)) peaks at w = ^ but the bracketed combination of functions vanishes, 
canceling the infrared divergence and leaving 



; m + #■ 



M{k^ 



|<:|<m 



m + # 



(5.9) 



so that from Eq. ( p.5| ), the functions and B,„ are given by 



MP) l + ^(l/m) 



(5.10) 



where for B,„ it is recognized that the infrared divergence of the integral is not suppressed by factors 
of l/m, and for A., we demand that m is the largest scale (the limit — oo is taken only at the end). 
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The spin-decomposed quark propagator is then 



The first component represents a heavy quark propagating forward in time (the second is the anti- 
quark propagating backwards in time) and explicitly agrees with the expression found in Ref. JI^]. 
Further, when considering the heavy quark limit of the Bethe-Salpeter equation, the Faddeev equa- 
tion or the quark four-point function, the constant and infrared divergent spatial integral occur- 
ring in B/, cancel explicitly when the quarks are in a color singlet configuration (and only for these 
configurations). The Bethe-Salpeter equation moreover furnishes the result that 

V{r)r^ J da>F{a>^){l-e"^ ') (5.12) 

(r is a length scale) so that the connection between the Coulomb kernel, F, and the quark-antiquark 
potential, V, is made explicit. 

The connection to the heavy quark limit is rather important, because it explains the role of the 
constant (arising from the charge constraint) and the infrared divergences. Within the leading 
order truncation presented here, the physical dynamics are contained within the static propagator 
dressing functions and their gap equations. The full propagators have pole positions that are de- 
pendent on and the infrared divergence. As the pole positions are shifted to infinity in the limit 
<^ —5. oo (i.e., when the total color charge is conserved and vanishing) and as the infrared integrals 
diverge, this simply reflects the fact that infinite energy is required to create isolated colored par- 
ticles from the colorless vacuum. However, for physical color singlets the divergent contributions 
cancel. It thus appears that the constant and the infrared divergence are merely constant shifts in 
the potential and which are not observable. 



6. Summary 

In summary, a leading order truncation to the Dyson-Schwinger equations of Coulomb gauge 
within the first order formalism has been presented. Because Coulomb gauge is incomplete, the 
temporal zero modes must be taken into account and it is seen that this results in a nonperturbative 
constraint on the total color charge. In the Coulomb gauge first order formalism, the ghosts cancel 
but the resulting action is nonlocal. To derive the Dyson-Schwinger equations, an Ansatz is thus 
made to replace the nonlocal Coulomb kernel with its expectation value. This introduces a set of 
four-point interaction vertices, which are dependent on the input Ansatz for the Coulomb kernel 
(F). Truncating the system to include only the tadpole diagrams involving this input Ansatz leads 
to a closed set of equations (their solution is discussed in Ref. Importantly, these equations 
reduce to the gap equations for the static gluon and quark propagators obtained from a quasi-particle 
approximation in the canonical Hamiltonian approach, Refs. [Q] and [Q], respectively. Furthermore, 
the known Coulomb gauge heavy quark limit [|l^] emerges. It is seen that the static propagator gap 
equations are not affected by the charge constraint or the infrared divergence of the input Ansatz, 
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F. This is in contrast to the full propagators, which are dependent. However, the connection 
to the heavy quark limit supplies an explanation: such unphysical singularities cancel for color 
singlet states, whereas the pole positions of colored propagators are shifted to infinity, reflecting 
that infinite energy is required for these to be created in isolation from the vacuum. 
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